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Показана возможность определения оптимального вида входного напряжения для электротехнических схем замещения перво- 
го и второго порядков барьерного разряда. Представлен алгоритм определения оптимальной формы воздействующего импуль- 
са напряжения, обуславливающего потребление минимальной энергии электротехническими элементами схемы замещения 
озонатора. 


В системах газо- и водоочистки большой инте- 
рес представляет использование импульсного 
барьерного разряда. Успешно осуществляемая в 
данном разряде реакция образования озона являет- 
ся одним из немногих плазмохимических процес- 
сов, реализованных в промышленном масштабе, 
причем мощности озонаторных установок непре- 
рывно растут. Этот тип разряда характеризуется, с 
одной стороны, сравнительно высокой средней 
энергией электронов 4.. .5 эВ и с другой - низкой 
температурой газа, которая близка к температуре 
электродов. При этом энергия, вложенная в раз- 
ряд, выделяется в короткоживущих малоинтенсив- 
ных искрах - микроразрядах. Сочетание всех этих 
условий делает барьерный разряд эффективным 
для осуществления реакций конденсации: и дей- 
ствительно, кроме уже упомянутого процесса обра- 
зования озона в нем можно проводить многие де- 
сятки органических и неорганических синтезов. 

При феноменологическом описании электри- 
ческих разрядов наибольшим успехом пользуется 
их описание как объектов электрической цепи 
[1-4]. В основе таких подходов лежит замена элек- 
трического разряда некой эквивалентной схемой 
замещения, состоящей из стандартных электротех- 
нических элементов. Особенности протекания фи- 
зических процессов определяются формой воздей- 
ствующего напряжения. В ряде случаев эффектив- 
нее оказывается использование импульсного на- 
пряжения питания. 

На рис. 1 приведена электротехническая схема 
замещения первого порядка для барьерного разря- 
да, предложенная авторами [2] . Суть данной рабо- 
ты заключается в том, что бы выбрать такую форму 
воздействующего напряжения, при котором вход- 
ной генератор выделял (затрачивал) бы наимень- 
шую энергию для достижения заданного значения 
напряжения на конденсаторе в фиксированный 
момент времени. Рассмотрим решение такой зада- 
чи для схемы, представленной на рис. 1. 

Для решения поставленной вариационной за- 
дачи [5, 6] необходимо определить переходную 
функцию по напряжению - И(1), связывающую 
входное напряжение и(1) с напряжением на кон- 
денсаторе и с (і). 

Переходная функция - реакция цепи на еди- 
ничное воздействие, поэтому, считая, что входное 
напряжение равняется единице, получаем: 


I 

/?(?) = 1-е (1) 




Рис. 1. Электротехническая схема замещения первого по- 
рядка разрядного промежутка озонатора 


Используя интеграл Дюамеля, можно записать 
напряжение на конденсаторе с учетом импульса 
воздействующего напряжения и полученной пере- 
ходной функции (1): 
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По условию задачи необходимо минимизиро- 
вать функционал энергии генератора: 
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при {= Т и дополнительном условии выполнения (2). 

Для вариации функционала необходимо выра- 
зить ток в цепи /(?), через входное напряжения и(1) 
и напряжение на конденсаторе и с (1). 


, (0 = » ( о-»ло 
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и, подставляя полученное соотношение (4) в выра- 
жение (3), с учетом(2) получаем: 
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Теперь можно осуществить условную миними- 
зацию функционала (5) с дополнительным услови- 
ем (2), которую формально можно записать так [5]: 
Ь = 1У + Л.-и с , 5Ь = 0, 


или с учётом фильтрующих свойств дельта-функции 


8(х-і) 


Г 0 при I Ф х; 
[оо при і = х 


и свойств функциональных производных [5]: 
3/((р(х)) _ д/((р(х)) 8(р(х) 8(р(х) _ 
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можно записать в развернутом виде: 
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После простых преобразований получаем, (по- 
яснения приводятся на рис. 2): 
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Рис. 2. 



Расстановка пределов интегрирования в кратном ин- 
теграле 


Таким образом, мы получили интегральное ура- 
внение Фредгольма второго рода, относительно 
неизвестной функции и(і) - входного напряжения: 
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Осуществляя алгебраизацию ур. (8), получаем си- 
стему линейных алгебраических уравнений относи- 
тельно искомой функции входного напряжения и(1) 
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- шаг дискретизации; 8, ре 0.../Ѵ; N - число точек 
дискретизации временного интервала. Запишем 
решение уравнения (9) в матричной форме, введя 
обозначения, Е={ 1 } - единичная матрица, 

(* = {&,*} - матрица проводимости, Ь={ — Я-/г А .} - ма- 
трица, образованная правой частью уравнения (9): 

К = С + — , -> Ѵ = К ь. 

м 

На рис. 3 приведён результат восстановления 
фронта импульса входного напряжения для мо- 
дельных параметров схемы К= 10 Ом; С=40 мкФ и 
длительности импульса Т=\ мкс. 

На рис. 4 представлены расчетные графики 
энергий для трёх разных форм импульсов входного 
напряжения (рис. 4 , а, б). Сплошной линией пока- 
зано входное напряжение - оптимальной расчёт- 
ной формы и соответствующая ему затраченная 
энергия генератора, пунктиром и штрих-пункти- 
ром изображены неоптимальные импульсы вход- 
ного напряжения и соответствующие им энергии. 
На графике видно, что энергия, выделяемая при 
оптимальной форме воздействующего напряже- 
ния, имеет меньшее значение по сравнению с 
энергиями, полученными в результате воздействия 
других импульсов. 



Рис. 3. Оптимальная форма импульса входного напряже- 
ния, полученная в результате расчёта 
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Рис. 4. Различные формы модельных импульсов входного 
напряжения и соответствующее им распределение 
энергии 


В ряде случаев необходимо учесть влияние про- 
водимости межэлектродной субстанции (1/К2 на 
рис. 5), например, в цепи второго порядка. Рассмо- 
трим алгоритм определения оптимальной формы 
воздействующего импульса, обуславливающей по- 
требление минимальной энергии электротехниче- 
скими элементами схемы замещения озонатора - 
цепи второго порядка на примере модельной зада- 
чи. Выберем в качестве величин цепи параметры: 
С1=510- 11 Ф, С2=410 10 Ф, Д1=15 Ом, Й2=15 Ом. 

Схема замещения барьерного разряда предста- 
влена на рис. 5, где С1 - ёмкость барьеров, ёмкость 
воздушной среды с вкрапленными каплями воды - 
С2, сопротивление схемы - К\, сопротивление 
межэлектродной субстанции - К1. 

Для решения поставленной вариационной зада- 
чи, как было показано выше, определим переход- 
ную функцию по напряжению - И(1), связывающую 
входное напряжение и(Г) с напряжением на конден- 
саторе и а (1) для схемы, изображённой на рис. 5. 

Считая, что входное напряжение равняется 
единице, получаем: 

Н(1)= А]в р ' 1 + А 2 е Р2 ‘ , (Ю) 
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А, + А= 0; 

А - А - 1 

А х р 1 +А 1 р 1 = — 1 — ; 1 ‘ {р х - р 2 )Я\- С2 

[ 11 22 К1-С2 


К2 



Рис. 5. Схема замещения барьерного разряда ~ цепь второ- 
го порядка 

Используя интеграл Дюамеля, можно записать 
напряжение на конденсаторе с учетом импульса 
воздействующего напряжения и полученной пере- 
ходной функции (10): 
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По условию задачи необходимо минимизиро- 
вать функционал энергии генератора: 


/ 
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( 12 ) 


при 1=Ти дополнительном условии выполнения (11). 

Для вариации функционала необходимо выра- 
зить ток в цепи і(і) через входное напряжение и(1) и 
напряжение на конденсаторе и а (і): 

і(1) = і, !2 (і) + і С2 (і) = С2* и “ (1) 
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Подставляя полученное соотношение в выра- 
жение (12), с учетом (11) получаем: 
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Теперь можно осуществить условную миними- 
зацию функционала (13) с дополнительным усло- 
вием (11), которую формально можно записать так: 

Ь = Ж + Л-и с , 81 = 0, 
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С учётом изложенного ранее, можно записать: 


8 _ 
8 и 


| Лм(/)| и(т)Н(1 -т)с/г 

о о 

Млг и(т)Н(\і-т\); 


К 


(14) 


8 _ 
8 и 


С і 


Я|и(т)Я(/-т)с/т = Я- Н(Т-1). (15) 


V о 


Выразим в виде суммы двух слагаемых: 
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Как было показано выше в (6-8) 
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Тогда результирующее выражение для миними- 
зируемого функционала с учётом (14-19) 
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Это интегральное уравнение Фредгольма вто- 
рого рода относительно переходной проводимости 
неизвестной функции входного напряжения и(Г): 
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Осуществляя алгебраизацию уравнения (20), 
получаем систему линейных алгебраических ура- 
внений относительно и(/) 
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- шаг дискретизации; $,р еО..Я; N - число точек 
дискретизации временного интервала. Запишем ре- 
шение уравнение (20) в матричной форме, введя 
обозначения, Е={ 1 } - единичная матрица, 

^ = {&} - матрица проводимости, Ь ={— Я-й*} - ма- 
трица, образованная правой частью уравнения (21): 

К = С + — Е, -> 11 = К ‘ Ь. 
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По аналогии с (14): 

К г г ^ 

С2р, | и(і)сіі^ и(т)<1т Я, 0 - г) 

ѵ о о 

т 

= С2р^Ат м(т)Я,(|/-т|); 
о 


8 _ 
8 и 


Таким образом, представленный авторами алго- 
ритм может быть использован для оптимизации ра- 
боты озонатора. Энергия, выделяемая генератором 
озонатора при оптимальной форме воздейству- 
ющего напряжения, имеет меньшее значение по 
( 18 ) сравнению с энергией, полученной в результате 
воздействия импульса неоптимальной формы. 
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